Chaste: an open source C++ library for computational physiology and biology

Chaste - Cancer, Heart And Soft Tissue Environment - is an open source C++ library for the computational simulation of mathematical models developed for physiology and biology. Code development has been driven by two initial applications: cardiac electrophysiology and cancer development. A large number of cardiac electrophysiology studies have been enabled and performed, including high performance computational investigations of defibrillation on realistic human cardiac geometries. New models for the initiation and growth of tumours have been developed. In particular, cell-based simulations have provided novel insight into the role of stem cells in the colorectal crypt. Chaste is constantly evolving and is now being applied to a far wider range of problems. The code provides modules for handling common scientific computing components, such as meshes and solvers for ordinary and partial differential equations (ODEs/PDEs). Re-use of these components avoids the need for researchers to "re-invent the wheel" with each new project, accelerating the rate of progress in new applications. Chaste is developed using industrially-derived techniques, in particular test-driven development, to ensure code quality, re-use and reliability. In this article we provide examples that illustrate the types of problems Chaste can be used to solve, which can be run on a desktop computer. We highlight some scientific studies that have used or are using Chaste, and the insights they have provided. The source code, both for specific releases and the development version, is available to download under an open source Berkeley Software Distribution (BSD) licence at http://www.cs.ox.ac.uk/chaste, together with details of a mailing list and links to documentation and tutorials

Larval Connectivity in an Effective Network of Marine Protected Areas

Acceptance of marine protected areas (MPAs) as fishery and conservation tools has been hampered by lack of direct evidence that MPAs successfully seed unprotected areas with larvae of targeted species. For the first time, we present direct evidence of large-scale population connectivity within an existing and effective network of MPAs. A new parentage analysis identified four parent-offspring pairs from a large, exploited population of the coral-reef fish Zebrasoma flavescens in Hawai'i, revealing larval dispersal distances ranging from 15 to 184 km. In two cases, successful dispersal was from an MPA to unprotected sites. Given high adult abundances, the documentation of any parent-offspring pairs demonstrates that ecologically-relevant larval connectivity between reefs is substantial. All offspring settled at sites to the north of where they were spawned. Satellite altimetry and oceanographic models from relevant time periods indicated a cyclonic eddy that created prevailing northward currents between sites where parents and offspring were found. These findings empirically demonstrate the effectiveness of MPAs as useful conservation and management tools and further highlight the importance of coupling oceanographic, genetic, and ecological data to predict, validate and quantify larval connectivity among marine populations

Transcriptomic identification of starfish neuropeptide precursors yields new insights into neuropeptide evolution

Published by the Royal Society under the terms of the Creative Commons Attribution License http://creativecommons.org/licenses/by/4.0/, which permits unrestricted use, provided the original author and source are credited.This work was supported by a PhD studentship funded by QMUL and awarded to D.C.S. and a Leverhulme Trust grant (RPG- 2013-351) awarded to M.R.E. Sequencing of the A. rubens neural transcriptome was funded by an EPSRC grant (EP/J501360/1

Stable Isotopic Evidence for Methane Seeps in Neoproterozoic Postglacial Cap Carbonates

The Earth's most severe glaciations are thought to have occurred about 600 million years ago, in the late Neoproterozoic era. A puzzling feature of glacial deposits from this interval is that they are overlain by 1–5-m-thick 'cap carbonates' (particulate deep-water marine carbonate rocks) associated with a prominent negative carbon isotope excursion. Cap carbonates have been controversially ascribed to the aftermath of almost complete shutdown of the ocean ecosystems for millions of years during such ice ages—the 'snowball Earth' hypothesis. Conversely, it has also been suggested that these carbonate rocks were the result of destabilization of methane hydrates during deglaciation and concomitant flooding of continental shelves and interior basins. The most compelling criticism of the latter 'methane hydrate' hypothesis has been the apparent lack of extreme isotopic variation in cap carbonates inferred locally to be associated with methane seeps. Here we report carbon isotopic and petrographic data from a Neoproterozoic postglacial cap carbonate in south China that provide direct evidence for methane-influenced processes during deglaciation. This evidence lends strong support to the hypothesis that methane hydrate destabilization contributed to the enigmatic cap carbonate deposition and strongly negative carbon isotopic anomalies following Neoproterozoic ice ages. This explanation requires less extreme environmental disturbance than that implied by the snowball Earth hypothesis

The catalytic subunit of the system L1 amino acid transporter (S<i>lc7a5</i>) facilitates nutrient signalling in mouse skeletal muscle

The System L1-type amino acid transporter mediates transport of large neutral amino acids (LNAA) in many mammalian cell-types. LNAA such as leucine are required for full activation of the mTOR-S6K signalling pathway promoting protein synthesis and cell growth. The SLC7A5 (LAT1) catalytic subunit of high-affinity System L1 functions as a glycoprotein-associated heterodimer with the multifunctional protein SLC3A2 (CD98). We generated a floxed Slc7a5 mouse strain which, when crossed with mice expressing Cre driven by a global promoter, produced Slc7a5 heterozygous knockout (Slc7a5+/-) animals with no overt phenotype, although homozygous global knockout of Slc7a5 was embryonically lethal. Muscle-specific (MCK Cre-mediated) Slc7a5 knockout (MS-Slc7a5-KO) mice were used to study the role of intracellular LNAA delivery by the SLC7A5 transporter for mTOR-S6K pathway activation in skeletal muscle. Activation of muscle mTOR-S6K (Thr389 phosphorylation) in vivo by intraperitoneal leucine injection was blunted in homozygous MS-Slc7a5-KO mice relative to wild-type animals. Dietary intake and growth rate were similar for MS-Slc7a5-KO mice and wild-type littermates fed for 10 weeks (to age 120 days) with diets containing 10%, 20% or 30% of protein. In MS-Slc7a5-KO mice, Leu and Ile concentrations in gastrocnemius muscle were reduced by ∼40% as dietary protein content was reduced from 30 to 10%. These changes were associated with >50% decrease in S6K Thr389 phosphorylation in muscles from MS-Slc7a5-KO mice, indicating reduced mTOR-S6K pathway activation, despite no significant differences in lean tissue mass between groups on the same diet. MS-Slc7a5-KO mice on 30% protein diet exhibited mild insulin resistance (e.g. reduced glucose clearance, larger gonadal adipose depots) relative to control animals. Thus, SLC7A5 modulates LNAA-dependent muscle mTOR-S6K signalling in mice, although it appears non-essential (or is sufficiently compensated by e.g. SLC7A8 (LAT2)) for maintenance of normal muscle mass

Characterization of New Substrates Targeted By Yersinia Tyrosine Phosphatase YopH

YopH is an exceptionally active tyrosine phosphatase that is essential for virulence of Yersinia pestis, the bacterium causing plague. YopH breaks down signal transduction mechanisms in immune cells and inhibits the immune response. Only a few substrates for YopH have been characterized so far, for instance p130Cas and Fyb, but in view of YopH potency and the great number of proteins involved in signalling pathways it is quite likely that more proteins are substrates of this phosphatase. In this respect, we show here YopH interaction with several proteins not shown before, such as Gab1, Gab2, p85, and Vav and analyse the domains of YopH involved in these interactions. Furthermore, we show that Gab1, Gab2 and Vav are not dephosphorylated by YopH, in contrast to Fyb, Lck, or p85, which are readily dephosphorylated by the phosphatase. These data suggests that YopH might exert its actions by interacting with adaptors involved in signal transduction pathways, what allows the phosphatase to reach and dephosphorylate its susbstrates

The orphan germinant receptor protein GerXAO (but not GerX3b) is essential for L-alanine induced germination in Clostridium botulinum Group II

Clostridium botulinum is an anaerobic spore forming bacterium that produces the potent botulinum neurotoxin that causes a severe and fatal neuro-paralytic disease of humans and animals (botulism). C. botulinum Group II is a psychrotrophic saccharolytic bacterium that forms spores of moderate heat resistance and is a particular hazard in minimally heated chilled foods. Spore germination is a fundamental process that allows the spore to transition to a vegetative cell and typically involves a germinant receptor (GR) that responds to environmental signals. Analysis of C. botulinum Group II genomes shows they contain a single GR cluster (gerX3b), and an additional single gerA subunit (gerXAO). Spores of C. botulinum Group II strain Eklund 17B germinated in response to the addition of L-alanine, but did not germinate following the addition of exogenous Ca2+-DPA. Insertional inactivation experiments in this strain unexpectedly revealed that the orphan GR GerXAO is essential for L-alanine stimulated germination. GerX3bA and GerX3bC affected the germination rate but were unable to induce germination in the absence of GerXAO. No role could be identified for GerX3bB. This is the first study to identify the functional germination receptor of C. botulinum Group II

Cluster analysis of protein array results via similarity of Gene Ontology annotation

BACKGROUND: With the advent of high-throughput proteomic experiments such as arrays of purified proteins comes the need to analyse sets of proteins as an ensemble, as opposed to the traditional one-protein-at-a-time approach. Although there are several publicly available tools that facilitate the analysis of protein sets, they do not display integrated results in an easily-interpreted image or do not allow the user to specify the proteins to be analysed. RESULTS: We developed a novel computational approach to analyse the annotation of sets of molecules. As proof of principle, we analysed two sets of proteins identified in published protein array screens. The distance between any two proteins was measured as the graph similarity between their Gene Ontology (GO) annotations. These distances were then clustered to highlight subsets of proteins sharing related GO annotation. In the first set of proteins found to bind small molecule inhibitors of rapamycin, we identified three subsets containing four or five proteins each that may help to elucidate how rapamycin affects cell growth whereas the original authors chose only one novel protein from the array results for further study. In a set of phosphoinositide-binding proteins, we identified subsets of proteins associated with different intracellular structures that were not highlighted by the analysis performed in the original publication. CONCLUSION: By determining the distances between annotations, our methodology reveals trends and enrichment of proteins of particular functions within high-throughput datasets at a higher sensitivity than perusal of end-point annotations. In an era of increasingly complex datasets, such tools will help in the formulation of new, testable hypotheses from high-throughput experimental data

Genome Wide Association Identifies PPFIA1 as a Candidate Gene for Acute Lung Injury Risk Following Major Trauma

Acute Lung Injury (ALI) is a syndrome with high associated mortality characterized by severe hypoxemia and pulmonary infiltrates in patients with critical illness. We conducted the first investigation to use the genome wide association (GWA) approach to identify putative risk variants for ALI. Genome wide genotyping was performed using the Illumina Human Quad 610 BeadChip. We performed a two-stage GWA study followed by a third stage of functional characterization. In the discovery phase (Phase 1), we compared 600 European American trauma-associated ALI cases with 2266 European American population-based controls. We carried forward the top 1% of single nucleotide polymorphisms (SNPs) at p<0.01 to a replication phase (Phase 2) comprised of a nested case-control design sample of 212 trauma-associated ALI cases and 283 at-risk trauma non-ALI controls from ongoing cohort studies. SNPs that replicated at the 0.05 level in Phase 2 were subject to functional validation (Phase 3) using expression quantitative trait loci (eQTL) analyses in stimulated B-lymphoblastoid cell lines (B-LCL) in family trios. 159 SNPs from the discovery phase replicated in Phase 2, including loci with prior evidence for a role in ALI pathogenesis. Functional evaluation of these replicated SNPs revealed rs471931 on 11q13.3 to exert a cis-regulatory effect on mRNA expression in the PPFIA1 gene (p = 0.0021). PPFIA1 encodes liprin alpha, a protein involved in cell adhesion, integrin expression, and cell-matrix interactions. This study supports the feasibility of future multi-center GWA investigations of ALI risk, and identifies PPFIA1 as a potential functional candidate ALI risk gene for future research

Novel variants in the PRDX6 Gene and the risk of Acute Lung Injury following major trauma

<p>Abstract</p> <p>Background</p> <p>Peroxiredoxin 6 (<it>PRDX6</it>) is involved in redox regulation of the cell and is thought to be protective against oxidant injury. Little is known about genetic variation within the PRDX6 gene and its association with acute lung injury (ALI). In this study we sequenced the <it>PRDX6 </it>gene to uncover common variants, and tested association with ALI following major trauma.</p> <p>Methods</p> <p>To examine the extent of variation in the <it>PRDX6 </it>gene, we performed direct sequencing of the 5' UTR, exons, introns and the 3' UTR in 25 African American cases and controls and 23 European American cases and controls (selected from a cohort study of major trauma), which uncovered 80 SNPs. <it>In silico </it>modeling was performed using Patrocles and Transcriptional Element Search System (TESS). Thirty seven novel and tagging SNPs were tested for association with ALI compared with ICU at-risk controls who did not develop ALI in a cohort study of 259 African American and 254 European American subjects that had been admitted to the ICU with major trauma.</p> <p>Results</p> <p>Resequencing of critically ill subjects demonstrated 43 novel SNPs not previously reported. Coding regions demonstrated no detectable variation, indicating conservation of the protein. Block haplotype analyses reveal that recombination rates within the gene seem low in both Caucasians and African Americans. Several novel SNPs appeared to have the potential for functional consequence using <it>in silico </it>modeling. Chi²analysis of ALI incidence and genotype showed no significant association between the SNPs in this study and ALI. Haplotype analysis did not reveal any association beyond single SNP analyses.</p> <p>Conclusions</p> <p>This study revealed novel SNPs within the <it>PRDX6 </it>gene and its 5' and 3' flanking regions via direct sequencing. There was no association found between these SNPs and ALI, possibly due to a low sample size, which was limited to detection of relative risks of 1.93 and above. Future studies may focus on the role of <it>PRDX6 </it>genetic variation in other diseases, where oxidative stress is suspected.</p